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Sideband holography can be used to get fields images ( Eo and E\) of a vibrating object for both 
the carrier (Eo) and the sideband (E i) frequency with respect to vibration. We propose here to 
record Eo and E\ sequentially, and to image the correlation E±Eq. We show that this correlation is 
insensitive the phase related to the object roughness and directly reflect the phase of the mechanical 
motion. The signal to noise can be improved by averaging the correlation over neighbor pixel. 

Experimental validation is made with vibrating cube of wood and with a clarinet reed. At 2 kHz, 
vibrations of amplitude down to 0.01 nm are detected. 
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There is a big demand for full field vibration measure¬ 
ments, in particular in industry. Different holographic 
techniques are able to image and analyze such vibrations. 
Double pulse jlj, Q or multi pulse holography Q records 
several holograms with time separation in the 1...1000 /is 
range, getting the instantaneous velocity from the phase 
difference. If the vibration frequency is not too high, one 
can also directly track the vibration of the object with 
fast CMOS cameras The analysis of the motion 

can be done by phase difference or by Fourier analysis in 
the time domain. For periodic vibration motions, mea¬ 
surements can be done with slow camera. Indeed, an har¬ 
monically vibrating object illuminated by a laser yields 
alternate dark and bright fringes [fi], that can be analyzed 
by time averaged holography [7]. Although the time av¬ 
eraged method gives a way to determine the amplitude 
of vibration [i[ quantitative measurement remains quite 
difficult for low and high vibration amplitudes. 

We have developed heterodyne holography 1, 0 , 
which is a variant of phase shifting holography, in which 
the frequency, phase and amplitude of both reference and 
signal beam are controlled by acousto optic modulators 
(AOM). Heterodyne holography is thus extremely versa¬ 
tile. By shifting the frequency of the local oscillator ujlo 
with respect to illumination uj o, if is f° r example possible 
to detect the holographic signal at a frequency co different 
than illumination ujq. This ability is extremely useful to 
analyze vibration, since heterodyne holography can de¬ 
tect selectively the signal that is scattered by object on 
vibration sideband of frequency = ujo + thuja , where 
uja is the vibration frequency and m and integer index. 

As was reported by Ueda et al, m the detection of the 
sideband m =1 is advantageous when the vibration ampli¬ 
tude is small. Nanometric vibration amplitude measure¬ 
ments were achieved with sideband digital holography 
on the m = 1 sideband 0, and comparison with single 
point laser interferometry has been made . Verrier et 


al. [HI have shown that one can simultaneously measure 
Eo and E\ by using a local oscillator with two frequency 
components. One can thus infer the mechanical phase 
of the vibration [15]. However, Bruno et al 0 shows 
that this simultaneous detection of Eo and E\ can lead 
to cross talk, and to a loss of detection sensitivity, which 
becomes annoying when the vibration amplitude is small. 

In this letter, we show that simultaneous detection of 
Eo and Ei is not necessary, and that equivalent or supe¬ 
rior performances can be obtained by detecting Eo and 
Ei sequentially. Indeed, the cross talk effects seen by 
Bruno et al. [H| disappear in that case. We also show 
that the random phase variations caused by the rough¬ 
ness of the object can be eliminated by calculating the 
correlation E\Eq. It is then possible to increase the sig¬ 
nal to noise ratio (SNR) by averaging correlation over 
neighboring pixels. The sequential measurement of Eo 
and E\, and the calculation of the correlation E\E^ 
make possible to image the vibration “full field”, and to 
measure quantitatively its amplitude and phase. Maxi¬ 
mum sensitivity is achieved by focusing the illumination 
in the studied point and by averaging the correlation in 
that region. Finally, we prove that the sensitivity is lim¬ 
ited by a sideband signal of one photo-electron per de¬ 
modulated image sequence. The device and method were 
validated experimentally by studying a cube of wood vi¬ 
brating at rz± 20 kHz, and a clarinet reed at ^ 2 kHz. 
Measurement sensitivity of 0.01 nm for a vibration at 2 
kHz, comparable to the sensitivity obtained by Bruno et 
al. [16| at 40 kHz, is demonstrated. 

Consider an object illuminated by a laser at frequency 
ujo that vibrate at frequency uja with an out of plane vi¬ 
bration amplitude z max . The out of plane coordinate is 
z(t) = z max sin (u)At)- The field scattered by the object 
is £ = E(t)e luJot + c.c., where c.c. is the complex conju¬ 
gate and E(t) the field complex amplitude. In reflection 
geometry, we have 

E(t) = E wo e Wsin{uAt+aTS ^ 

where E wo is the complex field without movement, 
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FIG. 1. Heterodyne holography setup applied to analyse vi¬ 
bration. L: main laser; AOM1, AOM2: acousto-optic modu¬ 
lators; M: mirror; BS: beam splitter; CCD: camera. 



FIG. 2. (a,b) Reconstructed images of a cube of wood vi¬ 
brating at uja / 27 t= 21.43 kHz. (a) £i£q correlation im¬ 
age: brightness is amplitude (i.e. |£i£q|), color is phase (i.e. 
aigE\E q). (b) 3D display of the phase arg£i£o- 


and $ is a complex quantity that describes the 
phase modulation. The phase of this modulation is 
arg<F, while the amplitude is |<F| = AnZmax/X. Be¬ 
cause of the Jacobi-Anger expansion, we have E(t) = 
Ewo £ m J m (|$|) e^ m ^ At+ar s^ where J m is the mth- 
order Bessel function of the first kind. The scattered 
field £ is then the sum of fields components £ m = 
E rn e luJrnt + c.c. of frequency uj rn = ojq + mu: a, where 
m is the sideband index (ra = 0 for the carrier) with 
£ m = E wo J m (|<f>|) e i mar s^. When the vibration am¬ 
plitude becomes small, the energy within sidebands 
decrease very rapidly with the sideband indexes m, and 
one has only to consider the carrier m = 0, and the first 
sideband mm z LI. We have thus: 

E 0 ($)=Ewo Jo(|$|) (1) 

E 1 (^)=E wo Ji(|$|) e iarg * 

Note that time averaged holography [7] that detects the 
carrier field £o is not efficient in detecting small ampli¬ 
tude vibration |<F|, because £o varies quadratically with 
|<F|. On the other hand, sideband holography that is able 
to detect selectively the sideband field E\ is much more 
sensitive, because E\ varies with linearly with |<F|. 

One must notice that the field scattered by the sam¬ 
ple without vibration E wo depends strongly on the x, y 
position. In a typical experiment, the sample rugosity is 
such that this field is a fully developed speckle. The re¬ 
constructed fields Eo and E\ are thus random in phase, 
from one pixel (x,y) to the next (x + l,y). One can¬ 
not thus simply extract the phase of the vibration from 
a measurement made on a single sideband. To remove 
the pixel to pixel random phase of E wo , we propose to 
record the hologram successively on the carrier m = 0 
and the sideband m — 1, to reconstruct the correspond¬ 
ing field image of the object Eo(x,y) and Ei(x,y), and 
to calculate and image the correlation £i£q, since this 
correlation do not involves E wo , but \E wo \ 2 , which is real 
and has no phase. Indeed, we have: 

( 2 ) 


For small vibration amplitude (|<F| 1), correlation sim¬ 

plifies to Ei Eq ~ \E wo \ 2 <$>/2. Correlation E\ £q is a 
powerful tool since gives directly the phase of mechani¬ 
cal motion arg<F. Nevertheless, problems can be encoun¬ 
tered when the signal |^7^ 0 1 2 scattered without vibration 
vanish. 

Figure [I] shows the holographic experimental setup 
used to measure successively £o and E\ in order to get 
E\Eq. The main laser L is a Sanyo DL-7140-201 diode 
laser (A=785 nm, 50 mW ). It is split into an illumination 
beam (frequency u;/, complex field £/), and in a LO beam 
(cj lo ,£lo). The illumination light scattered by the ob¬ 
ject interferes with the reference beam on the camera 
(Lumenera 2-2: 1616 x 1216 pixels of 4.4 x 4.4/im) whose 
frame rate is ujccd= 10 Hz. To simplify further Fourier 
transform calculations, the 1616x1216 measured matrix 
is truncated to 1024 x 1024. 

The illumination and LO beam frequencies ujlo and c01 
are tuned by using two acousto-optic modulators AOM1 
and AOM2 (Bragg cells), and we have ujlo = + 

uaomi and cj/ = cj l +^aom 2, where uaomi/2 —80 MHz 
are the frequencies of RF signals that drive the AOMs. 
The RF signals are tuned to have ujlo ~ ^1 — w ccd /4 
to get £ 0 , and to have ujlo — ^1 = + ^ccd /4 to get 

Ei. Successive sequences of Umax = 128 camera frames 
(i.e. Jo,/ 1 .../ 127 ) are recorded by tuning the RF signals 
first on the carrier (£ 0 ), then on the sideband (£ 1 ). The 
carrier and sideband complex hologram H are obtained 
from these sequences by 4 phase demodulation with n ma;E 
frames: 

Tl — Tl rn a x 1 

H(x,y)= j ni n(x,y) (3) 

n=0 

The fields images of the object Eo(x,y) and E\(x,y) are 
then reconstructed from H by the Schnars et. al Ga 
method that involves 1 Fourier transformation. The cor¬ 
relation £i£q is then calculated. 

Figure [2] (a) shows the reconstructed correlation im¬ 
ages of a cube of wood (2 cmx2 cm) vibrating at its res¬ 
onance frequency uja/2tt= 21.43 kHz. Brightness is the 


El E* = \E wo \ 2 Jid^l) Mm e iars * 
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correlation amplitude (i.e. |FiFg|) an d color the corre¬ 
lation phase (i.e. argFiFg). As seen, neighbor point 
have the same phase (same color). In order to get better 
SNR for the phase, the complex correlation signal E\E^ 
is averaged over neighbor x, y points by using a 2D Gaus¬ 
sian blur filter of radius 4 pixels. Figure [2] (b) displays 
the phase of the averaged correlation E\E^ in 3D. As 
seen, the opposite corners (upper left and bottom right 
for example) vibrate in phase, while the neighbor cor¬ 
ners (upper left and upper right for example) vibrate in 
phase opposition. Note that the cube is excited in one of 
its corner by a needle. This may explain why the opposite 
corners are not perfectly in phase in Fig. [2] (c). 

We can go further and use Eo and E\ to calculate the 
vibration complex amplitude <F. To increase the SNR 
let us average, over neighbor reconstructed pixels, the 
correlation E\E^ and the carrier intensity |Fo| 2 : 

\E 0 \ 2 = \E wo \ 2 J 0 2 ($) ~ \E wo \ 2 $/2 (4) 

By averaging we get: 

(EiE*) Xty = (1 /N pix ) Eli*', y’) E Si*', y’) (5) 

x',y' 

~ ($(x,y)/ 2) (1 /N pix ) \E wo {x’,y')\ 2 

x',y' 

(\E 0 \ 2 ) x ,y =* (1 /N pix ) Y \E wo {x’,y’)\ 2 

x',y' 

where ^2 X , , is the summation over the N p i x pixels of the 
averaging zone located around the point of coordinate 
x,y. The vibration amplitude <F is then 

<S>(x,y) = 2(E 1 E* 0 ) x J(\E 0 \ 2 ) x , y (6) 

We get here <f> that gives both the amplitude Zm ax = 
A|<F|/47r and the phase arg<f> of the mechanical motion. 
Note that it is also possible to get 4> by calculating the 
ratio Ei/Eo ~ <f> as done by Verrier et. al [l4|, but 
the ratio calculation is unstable for the points x, y of the 
object where E wo is close to zero. 

To evaluate the limits of sensibility of the correlation 
+ averaging method, we have calculated by Monte Carlo 
the detection limit of | <f> |, for an ideal holographic detec¬ 
tion that is only limited by shot noise. The calculation 
is similar to one made by Lesaffre et al. [18]. For each 
pixel (#, 7 /), each frame (n) and each sequence (m =0 or 
m=l), we have calculated the ideal camera signal I' n in 
the absence of shot noise. We have I' n = |^LO+j n ^o/il 2 ’ 
where the factor j n accounts for the phase shift of the lo¬ 
cal oscillator field Elo with respect to object field F 0 /i 
for frame n. To account for the roughness of the sample, 
E' 0 and E[ are taken proportional to a Gaussian speckle 
E' wo (x, y) that is uncorrelated from one pixel to the next, 
but remains the same for all frames n and all sequences 
m=0 or 1. To account for shot noise, we have added to I' n 
a Gaussian random noise s of variance of , where I' n 
is expressed in electron photo electron Units. The noise 
s is uncorrelated from one pixels x,y to another, from 



Total number of sideband photo electron: n ma xNpix(\Ei \ 2 ) 

FIG. 3. Ratio |(FlFo)|/(|Fo| 2 ) calculated by Monte Carlo 
by decreasing the sideband signal |Fl| 2 . Horizontal axis is 
the total sideband energy: nmaxN P i X {\E± | 2 ) in photo elec¬ 
tron Units. Simulation is made with Umax =400, N P i X — 50 2 , 
|Flo | 2 = 10 4 and |Fo| 2 = 10 2 photo electrons. 



FIG. 4. (a) Clarinet reed with illumination beam focused 

in xo,yo. (b) Sideband m — 1 reconstructed image of the 
vibrating reed. The display is made in arbitrary log scale for 
the field intensity \Ei(x,y)\ 2 . 


one frame n to another, and from one sequence m to an¬ 
other. By this way, we have obtained the Monte Carlo 
frame signals I n = I' n + s with whom we have performed 
the 4 phase demodulation with n ma;E frames of Eq. ©• 
We have then calculated the reconstructed signal Eq and 
Fi, the correlations and intensities EiEq and |Fo| 2 , and 
the means (EiEq) and (|F 0 | 2 ). We have then calculated 
the ratio (FiFo)/(|F 0 | 2 ) that gives 4> using Eq. (j 6]). 

Figure [3] gives the result of the Monte Carlo simulation 
for the ratio | (EiEq) |/(|F 0 | 2 ). Each point correspond 
to simulation made with a double sequence m=0 and 
1. The simulation is performed by decreasing the side¬ 
band averaged signal field intensity (|Fi| 2 )=1, 0.5, 0.25 
... photo electron per pixel and per frame. The other 
parameters of the simulation are n rnax = 400, N p i x = 50 2 , 
1 Elo | 2 =10 4 and (|Fo| 2 )=10 2 photo electrons. For each 
value of (|Fi| 2 ), 10 simulations are performed. As can be 
seen in Fig. [3] the ratio |(FiFo)|/(|F 0 | 2 ) decreases with 
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the sideband signal E\ proportionally with yj (|Fi| 2 ). 
When Ei becomes very small, the ratio reaches a noise 
floor related to shot noise. To simplify the present discus¬ 
sion, the results are displayed as a function of the total 
number of photos electrons on the sideband m = 1. The 
x-axis is thus n ma;E A/p^(|Fi | 2 4 ). As seen of Figj3]the noise 
floor is reached for n max 7Vp^(|Fi | 2 ) = 1. We have ver¬ 
ified by making simulation not displayed on Figj3] that 
this result do not depend on \Elo\ 2 and (|Fo| 2 ), pro¬ 
vided that \Elo\ 2 (|F 0 | 2 ) (|Fi| 2 ). The noise floor 

corresponds thus to a minimal ratio |(FiFo)|/(|Fo| 2 ) = 
1 /yj TimaxNpi x (| Eq | 2 ) (i.e. to 10 _4 ), and to a minimal de¬ 
tectable vibration amplitude <F = 2/yj Umax N p i x (|Fo| 2 ) 
(i.e. to 2 10 -4 ). 
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FIG. 5. Ratio \(EiEq)\/(\E o\ 2 ) as function of the reed 
excitation voltage converted in vibration amplitude Zmax in 
nm Units. Measurements: dark grey points, and theory 
Ji(|T|)/J 0 (|4>|) ~ |<f>|/2 : light grey curve. 


We have tested the ability of the correlation + av¬ 
eraging method to measure low vibration amplitude in 
an experiment made with a vibrating clarinet reed ex¬ 
cited at uja=2 kHz with a loudspeaker. In order to in¬ 


crease the sensitivity, the vibration amplitude is mea¬ 
sured on a single point (#o,2/o) ( see Tig0 (a)). The 
illumination has been focused on that point, and the 
calculations have been made with an averaging region 
centered on that point, whose size (radius 50 pixels for 
example) has been chosen to include the whole illumina¬ 
tion zone. We have reconstructed the field image of the 
reed at the carrier frequency (i.e. Eo(x,y)) (see Fig0] 
(b)), and verified on the holographic data that most of 
the energy |Fo | 2 is within the averaging zone xq,dq. Se¬ 
quences of Umax = 128 frames I n have been recorded for 
both carrier (Fo) and sideband (Fi), while the peak to 
peak voltage V pp of the loudspeaker sinusoidal signal has 
been decreased. We have then calculated H by Eq. ©, 
reconstructed the fields images of the reed Eq and Fi, 
and calculated the correlation FiFq, the intensity | Fo | 2 , 
and the means (FiFq) and (|F 0 | 2 ). We have then calcu¬ 
lated the ratio | (Fi Fq ) | / (| Fo | 2 ) - The latter is plotted on 
Fig|5] as a function of the loudspeaker voltage V pp that 
is proportional to vibration amplitude <F. The measured 
points follow |(FiFo)|/(|Fo| 2 ) ~ |4>|/2 oc V pp down to 
a vibration amplitude noise floor of |4>| ~ 7 10“ 5 that 
corresponds to 10 -2 nm i.e A/78000. Note that 

the noise floor measure here is about x20 lower than in 
previous holographic experiments [H, 141 and lower than 
the limit A/3500 predicted by Ueda LI]. Similar noise 
floor has been detected by holography by Bruno et al. 
[ill, but at much higher vibration frequency (40 kHz). 

In future work, it would be interesting to explore the 
limits of sensitivity of the correlation technique for low 
vibration amplitude. This can be done by using a laser 
with lower noise, by increasing the vibration frequency 
uja (and in all case by measuring the laser noise at uja)- 
One could also increase the illumination power, and the 
number of frame of the coherent detection (n ma:E > 128). 
A better control of the vibration uja versus camera ujccd 
frequencies, and or a proper choice of the demodulation 
equation could be also important, in order to avoid leak 
detection of Fo, when detection is tuned on E\. 
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